On the critical slowing down exponents of mode coupling theory 
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A method is provided to compute the parameter exponent A yielding the dynamic exponents of 
critical slowing down in mode coupling theory. It is independent from the dynamic approach and 
based on the formulation of an effective static field theory. Expressions of A in terms of third order 
coefficients of the action expansion or, equivalently, in term of six point cumulants are provided. 
Applications are reported to a number of mean-field models: with hard and soft variables and both 
fully-connected and dilute interactions. Comparisons with existing results for Potts glass model, 
ROM, hard and soft-spin Sherrington-Kirkpatrick and p-spin models are presented. 



In the framework of glassy mean-field (MF) models 
with quenched or built-in disorder, it is well known that 
the dynamics of models whose glassy phase is consistently 
described by a Replica Symmetry Breaking (RSB) solu- 
tion with a finite number of breaking's displays critical 
slowing down and a dynamic transition [1-4]. The equa- 
tions governing the relaxation dynamics down to the dy- 
namic critical temperature Td are those pertaining to the 
schematic mode-coupling-theory (MCT) developed in the 
context of supercooled liquids [5-9] . The "guide observ- 
able" is the correlation function C(t), i.e., the overlap 
between a given initial equilibrium configuration of the 
dynamics and the configuration at time t. In the dis- 
continuous one step RSB case near Td the relaxation of 
C(t) is a two step process in which the system spends a 
large amount of time, algebraically diverging as T — > Td, 
around a plateau value qea- Models with a discontinu- 
ous dynamic transition include, e.g., the spin-glass (SG) 
p-spin model with either spherical or Ising spins, the 
Potts glass model and the Random Orthogonal Model 
(ROM) . MCT predicts that two exponents control the 
whole dynamics. Exponent a governs the early ft regime 
as C(t) approaches the plateau value with a power law 
C(t) m qEA + c a /t a , while exponent b identifies the early 
a regime as C(t) = qEA — Cbt b , when the system starts re- 
laxing towards equilibrium [8]. An important prediction 
of MCT is the following relationship between the decay 
exponents and their mutual relationship to the so-called 
"exponent parameter" A: 

r 2 (i- a ) _ r 2 (i + 6) 
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In case of a continuous transition, there is no dynamic 
arrest and no b exponent is defined. Well known instances 
are, e.g., the paramagnet to full- RSB SG transition along 
the de Almeida Thouless (dAT) line in mean-field SG 
models, either fully connected [10, 11] or on random 
graphs [12], as well as the SG transition in Potts models 
with p < 4, both fully connected [13] and on the Bcthc 
lattice of any connectivity [14] , and in the p-spin spheri- 
cal model with large external magnetic field [3]. Even 
though Eq. (1) is usually assumed as the correct rela- 
tionship between exponents a and b, the exponent pa- 
rameter A, apart from schematic MCT cases, is simply 



considered a tunable parameter, without a specific con- 
nection to any physical observable. Here we show how to 
unveil this connection in full generality working within 
a "static-driven" effective theory of dynamics: we will 
put forward an independent formulation of A and apply 
it to some paradigmatic SG models. In general, the 
analytic treatment of the dynamics is more complicated 
than the statics and only a few models have been studied 
so far: the soft-spin Sherrington-Kirkpatrick (SK) model 
[10, 11], schematic MCT's [6], soft-spin p-spin and Potts 
glass models, for which the connection with MCT was 
first identified [2] . This prompted to consider the spheri- 
cal p-spin SG [3, 15] in all details as a MF structural glass, 
cf., e.g., Ref. [16], even in the off-equilibrium regime be- 
low Td [17]. In these cases dynamics is explicitly solved 
and A exactly computed. In particular, one finds that it 
is not universal and depends on model and external pa- 
rameters. On the other hand, its computation becomes 
difficult when we consider more complicated MF systems 
and finite-dimensional ones. 

Static definition of the exponent parameter A. 
Similarly to the static transition also the dynamic one 
can be located as the critical point of an appropriate 
replicated Gibbs free energy T [18-20]. This is a func- 
tion of the replicated dynamic variables (e.g. spins a), 
Legendre transform of the replicated free energy $ 



T[5Q ab ] = $[A Qb ] + A abSQ ab , 
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with respect to a conjugated field A, function of the over- 
lap matrix. We used the short-hand SQ a b = <^a<^b — q\ 
q = (o- a ab) being the Edwards Anderson parameter, that 
is the value of the overlap matrix elements for the replica 
symmetric (RS) solution. Average (. . .) is performed 
over the proper replicated ensemble [21] [46]. Expanded 
around the RS critical point T reads: 



T(SQ) 



with a,b,c,d = 1, 
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We have retained only two of all 



2 



cubic coefficients for they are the reievant ones at criti- 
cally [22] . The case n — is relevant for the continuous 
transition [23-25], the case n = 1 for the dynamic dis- 
continuous transition [18-20, 26]. Our main result is: 



A = 



W2 



(4) 



We note that this ratio also yields the breaking point x 
in the case of continuous transition from RS to full-RSB 
[13]. The above result can be obtained in the context of 
the supersymmetric formulation of the dynamics [27] that 
is a very convenient way of seeing the connection between 
equilibrium dynamics and the static replica treatment. 
Details of the derivation will be given elsewhere [28] . 

Eq. (4) holds in full generality above the upper crit- 
ical dimension. In general, we do not have an analytic 
expression of the Gibbs free energy, e.g., for MF models 
defined on finite-connectivity random graphs. However, 
r is defined as the Legendre transform of $ and, there- 
fore, its proper vertices can be associated to cumulants 
of the replicated order parameter. We, thus, face the 
problem of computing T from $ in presence of fields A ab 
coupled to l/Nj^i s i s i- The free energy $ needs to be 
computed at the third order, i.e., dealing with eight RS 
cumulants Wi,...,8 [22]. The coefficients w\ and w 2 of the 
Gibbs potential, cf. Eq. (3), however, are expressed as 



Wi 



r 3 uji 
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where r is the inverse of the SG susceptibility: 



(5) 
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i.e., the replicon, and u>\ and u>2 can be written in terms 
of six spin correlations (details elsewhere [28]): 



Wl = jy^ (siSj} c (sjS k }c(skSi)c 
ijk 
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where . . . means average over the disorder and (. . .) c 
denotes the connected thermal average. In case of dis- 
continuous transition it is implicit that different thermal 
averages in the above expression are all computed within 
the same state, selected by the initial condition. We un- 
derline that while vertices w\ and u>2 remain finite at T<j, 
the corresponding cu cumulants diverge as r~ 3 . 

Method validation and A computation. Applying 
Eq. (4) and computing w\ and W2 within the above 
mentioned static approach, we have tested our predic- 
tion on various models. In those few cases where A is 
exactly known from the dynamics [3, 11, 35] we verify 
that its formal expressions are identical. In more com- 
plicated systems the dynamic MCT phenomenology has 



been studied numerically and estimates of the exponents 
are available in the literature. These are instances in 
which Eq. 4 yields actual analytic prediction for A. 

Fully connected models. In the following, we consider 
a family of models with a Hamiltonian of the kind: 



H 



R(p) 



E JiiWi - E \/— E K,.,i P °n 



i<j 



p=l 



where <7j are N Ising spins, or soft /spherical ones. The 
2-body interaction matrix is constructed as J = T £0, 
where O is a random O(N) matrix chosen with the rota- 
tional invariant Haar measure, cf., e.g., Ref. [29] and £ 
is a diagonal matrix with elements independently chosen 
from a distribution p(e) [30]. In order to ensure the ex- 
istence of the thermodynamic limit, the support of p(e) 
must be finite and independent of N. The p-body inter- 
actions are i.i.d. Gaussian variables with zero mean 
and variance pl/iVP- 1 and = dP R{x) / dx? {x)\ x=0 for 
some real vauled function R(x). For the MF Ising SG 
[1, 4, 24], as well as for spherical SG's [15, 31, 32], the 
general form of the replicated free energy is: 

-n/3$ = extr QiA 5[Q,A] (9) 
S[Q,A] = lT T G(f3Q) + Lj2R(Q ab ) (10) 



--TrQA + lnTr w W[A;{<7}] 



W[A;{cr}] = cxp|i^A afc cr a cr fc | 



(11) 



where G : M n 



M nxn is a function in the space 



ofnxn matrices, formally defined through its power 
series around zero. Its form depends on the choice of the 
eigenvalue distribution p{e) of the £ matrix. 

Given this effective action S[Q,A], the saddle point 
equations in A and Q respectively read 



(12) 



Qab = (CaCh)w 

A ab = (3{G'(f3Q)} ab + l3 2 R'(Q ab ) 
that, in the RS Ansatz, become 
q = (m 2 (z)) z ; m(z) = (a) a (13) 
A = £[G' '(/3(1 + („-i)g))_G' (/3(1 -q))]+^Bf(q) 

where the weights over which (. . .) z and (...)„ are per- 
formed depend on a being Ising, spherical or soft. For 
the cases of interest here, weights are proportional to the 
following distributions: 



Ising 
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implying m(z) — tanh(z) for Ising and (1 — q)z for spheri- 
cal spins. To compute the values of Wi,2, cf. Eq. (3), one 
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2.05 


2.2 


2.5 
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5 
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A 


.5 


.556 


.652 


.719 


.743 


.746 


.743 


.739 


.736 


.733 


.731 


a 


.395 


.379 


.346 


.32 


.308 


.306 


.308 


.31 


.311 


.313 


.314 



TABLE I: Dynamic exponents in the Ising p spin model. 

has to expand Eq. (10) to third order around the saddle 
point value. Considering the second order term and im- 
posing that the Hessian determinant vanishes (criticality 
condition) we obtain 

(((a m f)\) z = [/3 2 G"(/3(1 - q)) + p 2 R"(q)] 1 (14) 

Using the above condition and expanding Eq. (12) to 
second order, we derive the general expressions of w's for 
fully-connected systems: 



3! 



/3 2 ~' 2 
3! 



G"'m-c i ))+A( q )({{a-mf)l) z (15) 
R"'(q) 



+ 2A{q){{{a-mY)i) 



A(q) = p 4 [G"(p(l-q))+R"(q)] i 



(16) 



holding both for n = and n = 1. They can be used to 
compute A in different cases. We now exemplify a few. 

For th e SK model , R{x) = x, G{x) = x 2 /2, p(e) = 
a/4 — e 1 1 (2-k), Sompolinsky's result for Ising spins along 
the dAT line is recovered [47], cf. Eq. 2.61 of Ref. [11]. 

For the ROM model [4, 30, 33], one has R{x) = 0, 



2G{x) = v a (x) - l + 2(a- 1) In 



v a (x) + 2x + 2a-l 



2a 



In 



v a (x) + 1 + 2x(2a: - 1) 



v a (x) = \j\ + Ax(2a - 1 + x) 

p(e) = aS(e — 1) + (1 — a)5(e + 1), and the transition is 
dynamic. MCT dynamics in the ergodic phase has been 
numerically studied in Ref. [34] for a = 13/32, where 
strong finite-size effects are observed and two different 
estimates for the exponent provided: b = .62, from the 
fit of the von Schwcidlcr law, while b — .75, from the fit of 
the equilibrium a relaxation time vs. temperature. Our 
"static-driven" computation yields b — .628 (A = .7077), 
allowing for a validation of the first numerical estimate. 

For Ising p-spin, G — and R(x) = x p /2, A values are 
in Tab. I. For p — > 2 we retrieve the result of Ref. [2]. 

For spherical models, G = 0, R(x) = ^ r x r /2, the 
exact analytic form A = k"(qd){l — qd) 3 /2 is retrieved 
[3, 35, 36], equivalent to the schematic MCT prediction 
for very long times at criticality, i.e., for C(t) = qd [7, 9], 
A(C(t)) being the MCT equations memory kernel. 

Yet another example is the Potts glass: 



ft = - X] Ji i ( p ~ l ) 



(17) 



The model has a discontinous glass transition for p > 4 
[13]. Brangian et al. studied the Potts glass with p — 10 



by means of Monte- Carlo numerical simulations [37] . Ap- 
proaching the dynamic transition, finite-size effects turn 
out to be large, implying that the plateau is almost in- 
visible also for very large sizes: this makes the numerical 
estimation of exponents very difficult. Their interpola- 
tion yields a = .33 ± .04. For p = 10, from the expan- 
sion of T[SQ] around q = qd, we obtain the exact values 
A = .8053 and a = .2759 [38], compatible with, though 
not extremely near to, the numerical estimate. 

Models on diluted random graphs. To study glassy 
models on random graphs we set up an apart technical 
method to analytically compute dynamic exponents. To 
frame the results, we first recall that there is no single 
instance of these class of models whose dynamics has 
been solved explicitly. Eq. (4) would allow to bypass 
the problem provided one had the replicated action T. 
Unfortunately, the action is hard to compute in diluted 
models, thus preventing w's derivation. We can, however, 
explicitly compute the related cumulants Wi ; 2 at critical- 
ity, cf. Eqs. (7)-(8). Details of the derivation will be 
given elsewhere [38], while here we present a validation 
of the results for the Viana-Bray model [39] in a field, 
displaying a continuous transition along the dAT line. 

Let us consider a Bethe lattice with connectivity c = 4, 
in a field h = .7 at the corresponding critical temper- 
ature T c (h) = 0.73536(1) [40]. The cumulants ratio 
yields A = 0.461. Our analityc prediction for the j3 de- 
cay exponent is, thus, a t h = .406. Wc, then, numer- 
ically study the system by means of both equilibrium 
and off-equilibrium Monte Carlo simulations. At equi- 
librium sizes 2 8 < N < 2 13 have been probed. Starting 
from an equilibrated configuration we measure the decay 
of C(t) = J2i s i(®) s i(t)/N with time. For each sample 
6 replicas are simulated and 16 uncorrelated measure- 
ments taken for each replica. Number of samples: from 
1280 (N — 2 13 ) to 5120 (N = 2 6 ). Off-equilibrium, the 
two time C(t,t w ) is measured. For very long t w , C(t,t w ) 
tends to the equilibrium C(t). For each sample we take 
one measure (at the longest simulated t w ). The number 
of samples goes from 1280 (L = 2 20 ) to 51200 (L = 2 6 ). 
Due to finite size effects, C(t) (at and off- equilibrium) 
displays a power law behavior only for times smaller than 
a time scale t*(N) that diverges with N. This makes the 
estimation of the exponent hard. Therefore, rather than 
C(t), we probe its thermal and disorder fluctuations, 



X4(t)=N (C(tr) - (C(t)) 



(18) 



remaning finite for large N at finite t. Using finite-time 
scaling arguments it can be argued that (i) \4 (*) diverges 
as t a at large times on a scale smaller that t*(N), (ii) the 
critical region diverges as t*(N) = N 1 ^ 3 ^, (iii) on times 
scales larger than the critical region the fluctuations scale 
as iV 1 / 3 [41]. As a consequence, if a has the correct value, 
the rescaled dynamic critical Xi(t)/ Nl/3 vs - t/A^ 1/(3a) , 
should be size independent. We plot it in Fig. 1 with 
a — a t h = 0.406. Collapse appears excellent both at and 
off-equilibrium. 
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FIG. 1: Rescaled X^if) f° r the Viana-Bray model with c = 4, 
h = .7 and T = .73536(1), along the dAT line with a = a th = 
0.406. No fitting interpolation is performed. 



Concluding, we have introduced a "static-driven" 
method to obtain, by means of a replica field theory, 
the dynamic exponents of the critical slowing down. The 
method allows to determine the MCT parameter expo- 
nent A as the ratio of coefficients of third order terms of 
the Gibbs free energy action expanded around the criti- 
cal point, cf. e.g., Eqs. (15-16). Equivalently, A is shown 
to be equal to the ratio between six points cumulants of 
a theory whose action is the Legendre transformed of the 
Gibbs free energy, cf. Eqs. (7-8). Indeed, the dynami- 
cal exponents can be associated to the ratio between two 
physical observables computed within a static framework. 
We verified the method's prediction in various MF mod- 
els, both on fully-connected and diluted graphs, success- 
fully comparing with previous analytical and numerical 
results. The method can be exported to any glass mod- 
els whose Gibbs action is computable or whose six point 
cumulants can be estimated. 
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